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Abstract. We study the generation of a stochastic gravitational wave (GW) background 
produced by a population of neutron stars (NSs) which go over a hadron-quark phase transition 
in its inner shells. We obtain, for example, that the NS phase transition, in cold dark 
O matter scenarios, could generate a stochastic GW background with a maximum amplitude 

of /ibg ~ 10 24 , in the frequency band ~ 20 — 2000 Hz for stars forming at redshifts of up to 
5^ z ~ 20. We study the possibility of detection of this isotropic GW background by correlating 

signals of a pair of 'advanced' LIGO observatories. 
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1. Introduction 

It is widely accepted that neutron stars (NSs) are born with fast rotating angular velocities. 
Due to magnetic torques, however, the NS periods could well be spun down. This spin-down 
causes a reduction in the centrifuge force and, consequently, the central energy density of the 
' NSs increases. Those stars, born with densities close to that of the quark deconfinement, may 

undergo a phase transition forming a strange quark matter core. It is worth stressing, however, 
that it is not a common sense that such a transition really takes place. 

As a consequence of such a putative phase transition, the stars could suffer a collapse which 
could excite mechanical oscillations. The great amount of energy generated in this process, 
AE ~ 10 53 ergs (~ 0.1 M & c 2 ) pQ, could be dissipated, at least partially, in the form of 
gravitational waves (GWs), which are studied in the present work. 

We adopted in the present paper the history of star formation derived by Springel & Hernquist 
[2], who employed hydrodynamic simulations of structure formation in a A cold dark matter 
(ACDM) cosmology. These authors study the history of cosmic star formation from the "dark 
ages", at redshift z ~ 20, to the present. 

Besides the reliable history of star formation by Springel & Hernquist, we consider the role of 
the parameters aj=i 2, which gives the fraction of the progenitor mass which forms the remnant 
NSs. We consider that the remnant mass is given as a function of the progenitor mass, namely 
M r = aim + QJ2- Later on we justify why M r is written in this way. Recall that a given initial 
mass function (IMF) refers to the distribution function of the stellar progenitor mass, and also 
to the masses of the remnant compact objects left as a result of the stellar evolution. 

In the present study we have adopted a stellar generation with a Salpeter IMF, which is 
consistent with Springel & Hernquist, since they show that population II stars could have been 
formed at high redshift too. We then discuss briefly what conclusions would be drawn whether 



(or not) the stochastic background studied here is detected by GW observatories such as LIGO 
and VIRGO. 

The paper is organized as follows. Section 2 deals with the NS EOSs adopted in our study; 
in Section 3 we consider the IMF and the NS masses; Section 4 deals with the GW production; 
in Section 5 we present the numerical results and discussions; Section 6 the detectability of the 
background of GWs are considered and finally in Section 7 we present the conclusions. 

2. Neutron stars and the nuclear matter equation of state 

We consider in the present study a particular model developed by Taurines et al [3] for a new 
class of parameterized field-theoretical model described by a Lagrangian density, where the whole 
baryon octet is coupled to scalar and vector fields through parameterized coupling constants. 
The free lepton fields contributes to the electrical equilibrium in the NS matter. 

Using a parametrization for the baryon-meson coupling constants, this model describes a wide 
range of NS parameters, such as a maximum mass ranging from (very low values) M = 0.66M© 
up to M = 2.77 'Mq, and the corresponding radii that vary in the range of 8 < R < 13 km. Each 
pair in the mass-radius relation is associated with a different parametrization of the EOS (see 
|3j for further details). Of course, each value of such parametrization represents different values 
of nuclear matter properties. 

Another scenario to be discussed concerns the strange quark star. The model used in this 
analysis is the MIT bag model [4J. The effects of pressure difference in the interior and exterior 
regions of the bag are summarized in the bag constant. 

Once the transition to quark-gluon matter occurs, the weak interaction processes for the 
quarks u, d and s 

u + s — > d + u (1) 

and 

d + u^ u + s (2) 

will take place, and a rapid transition occurs with a consequent gravitational micro-collapse. 
As a result, a huge amount of energy is dissipated, some in the form of GWs produced by 
quasi-normal modes excitation. 

3. Initial mass function and the neutron star masses 

The calculation of the GW background from NS phase transition requires the knowledge of the 
distribution function of stellar masses, the so called stellar initial mass function (IMF), <f)(m). 
Here the Salpeter IMF is adopted, namely 

<j)(m) = Am-( 1+X \ (3) 

where A is the normalization constant and x = 1.30 (our fiducial value). The normalization of 
the IMF is obtained through the relation 

rrriu 

I m(f)(m)dm = 1, (4) 

J JTtl 

where we consider m\ = 0.1 Mq and m u = 125 M . For further details we refer the reader to, 
e.g., [5]. 

It is worth mentioning that concerning the star formation at high redshift, the IMF could be 
biased toward high-mass stars, when compared to the solar neighborhood IMF, as a result of 
the absence of metals [H [7] . 

On the other hand, in the study by Springel and Hernquist, besides a high- mass star 
formation, it is shown that the population II stars, whose IMF could well be of the Salpeter's 
type, could start forming around redshift 20 or higher. 



Table 1. The values of oti=\ 2 which determines the fraction of the NS remant. 



a x 1/32 1/50 1/17 

a 2 (M Q ) 3/4 1.1 0.53 



In the present study we consider this population II studied by these authors. Then, for the 
standard IMF, the mass fraction of NSs produced as remnants of the stellar evolution is 

rrriu 

/ns= / M t </>(m)dm, (5) 

'-' m min 

where m m j n is the minimum stellar mass capable of producing a NS at the end of its life, and M r 
is the mass of the remnant NS. Stellar evolution calculations show that the minimal progenitor 
mass to form NSs is m m i n = 8M©, while the maximum progenitor mass is m max = 25 — 40M Q 
(see, e.g., 0). 

We do not intend to discuss in the present paper stellar evolution scenarios related to the 
formation of NSs. Many works can be found in the literature concerning the study of stellar 
evolution, supernova and NS formation [9]. However, we need to formulate a relation between 
the progenitor star mass m and the remnant NS mass M r . For the remnant, M r , we take 



M r = ai m + a 2 , 

where a\ and a 2 are constants; the first one is dimensionless and the second one is given in 
solar masses. As the results for stellar evolution are not yet fully determined, we studied this 
parametrization in three different scenarios for the values of a± and a 2 which would represent 
the distribution of NS mass as function of the progenitor mass (see Table 1 and the following 
sections). For example, considering NSs formed from progenitors with m = 20M Q we would 
have remnants with masses M r = 1.375M0, 1.5M Q and 2.7 Mq. These parameters describe 
sharper or softer distribution of NS masses around the standard value of mass M = 1.4M . 
Note that we are considering in the end that M r is independent of the redshift. 

With these considerations at hand, the mass fraction of NSs reads up to /ns = 10~ 2 for 
x = 1.30, while the fraction of NSs that undergoes a phase transition (/^g) can drop down to 
values <C 1 (see, e.g., [TO]). 



4. Gravitational wave production 

The dimensionless amplitude of the Gravitational Wave Background from Neutron Star Phase 
Transition can be given by 



h 2 



— I h 2 NS dR, (6) 

'obs J 



Vobs 

(see HBEl]). 

The micro-collapse produces GWs at frequency v of the NS f-mode, and dimensionless 
amplitude given by [13] 

E X 1 / 2 {2kHz\ 1/2 flMpc\ 



where E is the available pulsation energy, and d-^ is the luminosity distance to the source. 

Recall that the energy available to excite the pulsating modes is directly related to the EOS 
adopted. 



For the differential rate of NS micro-collapse we have 



di?NS = P*{z) — f^ s cf)(m)dmdz 



(8) 



where p*(z) is the star formation rate (SFR) density (in M yr _1 Mpc -3 ), dV is the comoving 
volume element, and /j^g, as already mentioned, is the fraction of NSs which may undergo a 
phase transition forming a strange quark matter core. 

For the SFR density, we adopt the one derived by Springel and Hernquist (see, [2] for details), 
namely 



where A = 3/5, /3 = 14/15, 

z m — 5.4 marks a break redshift, and Pm = O.^Moyr^Mpc" 3 
fixes the overall normalization. 

It is worth mentioning that these authors employed hydrodynamic simulations of structure 
formation in a ACDM cosmology with the following parameters: = 0.3, = 0.7, Hubble 
constant Hq = 100 h km s _1 Mpc -1 with h = 0.7, Qb = 0.04, and a scale-invariant primordial 
power spectrum with index n = 1, normalized to the abundance of rich galaxy clusters at present 
day (ct 8 = 0.9). 

Another relevant physical quantity associated with the GW background, produced by the 
first stars, is the closure energy density per logarithmic frequency span, which is given by 



In the next section we present the numerical results and discussions, which come mainly from 
the equation for /ibg- 

5. NUMERICAL RESULTS AND DISCUSSIONS 

In order to cover a wide number of parameters we present, in Table 2, the models considered in 
our study. In the first column appears the name of the model; in the second and third columns 
we present different combinations of the parameters a\ and 02, which imply in different ways 
to calculate the NS remnant mass for a given IMF; in the fourth column the mass range of the 
progenitor star; in the fifth column the NS remnant mass; in the sixth column the NS redshift 
formation, and finally in the seventh column the observed frequency, which are obtained via the 
empirical formulae given in reference }14j . 

These parameters are directly related to the number of NSs that goes over the phase 
transition. According to these NS models, only a fraction of them develops a core composed 
by deconfined quarks and gluons. When we set the mass range of the progenitor mass, we are 
also setting the mass of the remnant object. Stars with higher mass reach central densities that 
are high enough to develop the quark core, while smaller stars do not. The available energy for 
GW emission generated in the collapse process is also implicit in the mass of the remnant star. 
The more massive is the star, greater is its core and greater is the energy difference between the 
neutron and the hybrid star. We have found that stars with baryonic mass Mj, = 2M & develops 
a core as large as R = 7km and can generate as much as 3 x 10 53 erg of energy, while stars with 
baryonic mass M\, = 1M develops a small core with R = lkm and has ten times less energy 
available from the transition. As we do not know for sure the amount of this energy that is 
driven in each mode and in order to simplify our calculations we have adopted the standard 
medium value of O.OIM© for the released energy. 

In figure Q] we compare the background of GWs generated by NSs, which undergo phase 
transition, with that by the BH formation, which undergo quasi-normal mode instability (see 




0) 



]__dpcw_ 

p c dlogi^ obs ' 



(10) 



Table 2. Models description: mass function constants, ai2> progenitor, m, and remnant, M, 
mass ranges, redshift, z, and the observed frequency, v. 



Model 


OL\ 


02 


Am(M ) 


AM(M ) 


Az 


Au(Hz) 


A 


1/32 


3/4 


8.0-40.00 


1.00-2.00 


0-20 


20-2000 


B 


1/32 


3/4 


14.4-16.32 


1.20-1.26 


0-20 


31-1666 


C 


1/32 


3/4 


14.4-16.32 


1.20-1.26 


0-20 


75-1666 


D 


1/32 


3/4 


20.8-22.72 


1.40-1.46 


0-20 


27-1428 


E 


1/50 


1.1 


8.0-40.00 


1.26-1.90 


0-20 


21-1587 


F 


1/50 


1.1 


15.0-18.00 


1.40-1.46 


0-20 


27-1428 


G 


1/17 


0.53 


8.0-25.00 


1.00-2.00 


0-20 


20-2000 



|15j for details). In particular, for the NSs we consider two situations. We refer the reader to 
Table 2 to see the parameters adopted in these calculations. 

le-22 p , , 1 z 




v(Hz) 

Figure 1. The dimensionless amplitude for the GW background of stellar black hole formation 
(solid line) , the whole ensemble of NSs undergoing phase transition (dotted line) and 1 % of the 
NSs doing so (dashed line), for < z < 20. 

We consider in figure Q] two versions of the model B: one with its original values, namely, 
/ns = 0-01 and the other one, just as an example, with f^ s = 1. Obviously, even the most 
favorable values for nuclear and subnuclear matter coupling constants could not make even 
tend to 1, but this model is included just as a matter of comparison, as it would represent an 
upper limit. In both cases we consider E=0.01, i.e., only a tenth of the energy generated in 
phase transition goes to the pulsating mode. 

Note that the background of GWs generated by the NSs would have an amplitude greater 
than that generated by the BHs only if a considerable fraction of the NSs formed undergo phase 
transition. 

We compare how different laws to calculate the NSs mass, the remnant mass, for a Salpeter 
IMF, modify the spectrum of the background of GWs generated. We obtain that there is no 
significantly difference in the background for the three cases considered. 

It is worth mentioning that a detailed discussion of the models presented in Table 2, among 
other issues, will be considered in another publication to appear elsewhere. 



Table 3. For the models of Table 2 with different values of ffjg, we present the S/N for pairs 
of initial, enhanced and advanced LIGO observatories for one year of observation. We consider 
in the calculations E = O.OIMqc 2 . 





f pt 
Jns 




S/N 




Model 


Initial LIGO 


Enhanced LIGO 


Advanced LIGO 


A 


1.0 


3.6 x 10~ a 


1.8 x KT 1 


4.3 x lO" 1 


B 


0.01 


6.5 x 10~ 6 


2.6 x 10" 4 


5.0 x 10~ 4 


C 


0.01 


6.5 x 10" 6 


2.6 x 10~ 4 


5.0 x 10~ 4 


D 


0.01 


7.5 x 10~ 6 


3.1 x 10 -4 


6.5 x 10~ 4 


E 


1.0 


2.8 x 10 -3 


1.4 x io _1 


3.3 x lO" 1 


F 


1.0 


1.7 x 10~ 5 


7.0 x 1(T 4 


1.5 x 10~ 3 


G 


1.0 


2.7 x 10~ 3 


1.4 x 10" 1 


3.3 x 1Q- 1 



6. DETECTABILITY OF THE BACKGROUND OF GRAVITATIONAL WAVES 

The background predicted in the present study cannot be detected by single interferometric 
detectors, such as VIRGO and LIGO (even by advanced ones). However, it is possible to 
correlate the signal of two or more detectors to detect the background that we propose exists. 

To assess the detectability of a GW signal, one must evaluate the signal-to-noise ratio (S/N), 
which for a pair of interferometers is given by (see, e.g., [161 fT7] ) 



(S/N) 2 




dv 



(11) 



where sj® is the spectral noise density, T is the integration time and is the overlap reduction 
function, which depends on the relative positions and orientations of the two interferometers. 
For the y(v) function we refer the reader to |16j . who was the first to calculate a closed form for 
the LIGO observatories. 

Here we consider, in particular, the LIGO interferometers. Their spectral noise densities 
have been taken from [18] - who in turn obtained them from Thorne, by means of private 
communication. It is worth mentioning that although a ten years old reference has been cited, 
it presents analytical equations for the spectral noise densities, which are in complete agreement 
with many recent papers involving detectability of GWs by the LIGOs. 

In Table 2, as already mentioned, we present the models considered in our study. We show 
in Table 3 the S/N for the models of Table 2 for the three different LIGO generations and for 
different values of fjyg. 

As shown in Table 3, the signal-to-noise ratio for all models studied is lower than one, even 
for an advanced LIGO. Therefore, contrary to the claim of [19] such a putative GW background 
would hardly be detected. 

Note that the signal-to-noise ratio, for given IMF, Oj, and integration time, depends on 
and E as follows 

(S/N) oc (12) 

and it also depends on the SFR density in a more complicated way, namely, through an integral 
involving the redshift z. The higher the star formation rate, the higher the signal-to noise ratio 
will be. 

Just a matter of comparision, even considering models with f^ s = 1, a signal-to- noise ratio 
significantly greater than one, for the advanced LIGO, would be possible either the SFR density 



would be much greater than that by Springel & Hernquist, or the energy generated in the phase 
transition were almost completely channeled to excite the f- mode. Obviously, an optimist 
combination of the SFR density and the energy channeled to the f-mode would also render the 
same. 



7. CONCLUSIONS 

We present here a study concerning the generation of GWs produced from a cosmological 
population of NSs. These stars may undergo a phase transition if born close to the transition 
density, suffering a micro-collapse and exciting quasi-normal modes. 

We show that a detectable background is possible only if the SFR density is much greater 
than that predicted by Springel & Hernquist or if the energy generated in the phase transition 
is almost completely channeled to excite the f- mode. Obviously, a too optimistic combination 
of these two possibilities could do the same. 
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